Magnetic Structure of Hydrogen Induced Defects on Graphene 



J. O. Sofo, 1 Gonzalo Usaj, 2 P. S. Cornaglia, 2 A. M. Suarez, 1 A. D. Hernandez-Nieves, 2 and C. A. Balseiro 2 

'Physics Department, Pennsylvania State University, University Park, PA 16802, USA 
2 Centro Atomico Bariloche and Instituto Balseiro, CNEA, 8400 Bariloche, and CONICET, Argentina 

(Dated: submitted, December 201 1) 

Using density functional theory (DFT), Hartree-Fock, exact diagonalization, and numerical renormalization 
group methods we study the electronic structure of diluted hydrogen atoms chemisorbed on graphene. A com- 
parison between DFT and Hartree-Fock calculations allows us to identify the main characteristics of the mag- 
netic structure of the defect. We use this information to formulate an Anderson-Hubbard model that captures the 
main physical ingredients of the system, while still allowing a rigorous treatment of the electronic correlations. 
We find that the large hydrogen-carbon hybridization puts the structure of the defect half-way between the one 
corresponding to an adatom weakly coupled to pristine graphene and a carbon vacancy. The impurity's magnetic 
moment leaks into the graphene layer where the electronic correlations on the C atoms play an important role 
in stabilizing the magnetic solution. Finally, we discuss the implications for the Kondo effect. 

PACS numbers: 73.22.Pr,81.05.ue,73.20.-r,73.20.Hb 



I. INTRODUCTION 

The unusual electronic and mechanical properties of 
graphene 1-6 triggered an intense research activity after the 
first isolation of individual graphene sheets. 1,7 ' 8 The anoma- 
lous charge transport properties of pure graphene are a con- 
sequence of its particular sub-lattice symmetry that leads to 
a band structure with a point-like Fermi surface and a linear 
dispersion close to the Fermi energy, 9 These properties are 
responsible for the peculiar low energy excitations, that cor- 
respond to massless-Dirac fermions, and for the long electron 
mean-free path. 

The possibility of doping graphene with electrons or holes 
using external gates, which does not lead to a significant 
loss of mobility, allows the design of graphene-based de- 
vices. Another approach, that is also viewed as a promis- 
ing way to engineer the band structure and modify its elec- 
tronic properties, is the controlled adsorption of different 
types of atoms or molecules on grapheneJ^— In particu- 
lar, there has been an increasing interest to incorporate mag- 
netic effects in graphene to use graphene-based devices in 
spintronics. 14 ' 19 ! 26- — There are already some advances in this 
direction, like the injection of spin polarized currents from fer- 
romagnetic electrodes i 26 ' 31 ' 32 An alternative route is the gen- 
eration of magnetic defects using, for example, adatoms, va- 
cancies or edges. Hydrogen (H) impurities and carbon (C) 
vacancies are among the simplest and most studied point-like 
magnetic defects^— 

The problem of H atoms chemisorbed on graphene has 
been studied by several groups i 11 ' 16 ' 17 ^ 3 ' 30 ' 39,40 It is known 
that H impurities are adsorbed on top of a C atom, form- 
ing a covalent bond and locally distorting the honeycomb lat- 
tice. There is an increasing consensus that H adatoms acquire 
a magnetic moment, therefore behaving as magnetic impu- 
rities. A simple picture used to describe H on graphene is 
based on the Anderson model where the localized Is H orbital 
is hybridized with the graphene extended states. The intra- 
atomic Coulomb repulsion at the ls-orbital, together with the 
graphene pseudo-gap favors the existence of a local moment 
at the H impurity 



The case of C vacancies is quite different. In the ideal case 
where the lattice remains undistorted after the removal of a 
C atom, 27 there is a reduction of the coordination number of 
the three neighboring atoms that generates a resonant state 
at the pseudo-gap.— The existence of this resonance together 
with the intra-atomic Coulomb repulsion at the C orbitals nat- 
urally leads to the formation of a magnetic moment. This 
simple picture is consistent with the observation of magnetic, 
Kondo-like correlations in irradiated samples 42 and enlightens 
the role of the on-site Coulomb interaction in graphene (from 
hereon denoted as U~c)- The parameter Uc has been estimated 
to be of the order of a few electronvolts. 43 

A realistic description of impurities on graphene should 
then include the effects of the electronic correlations in both 
the impurity and the C orbitals by means of an Anderson- 
Hubbard like model. Such a model can extrapolate between 
two different limits: (z) for a small value of the hybridization 
between the impurity and C orbitals, a magnetic moment that 
is mainly localized at the impurity orbital can be stable; and 
(z'z) for a large hybridization, bonding and anti-bonding states 
are shifted away of the low energy region. The effect of these 
shifts is to effectively remove the p z orbital of the hybridized 
carbon atom from the Fermi energy generating a vacancy-like 
defect with a magnetic moment localized at the C atoms. 

As we will discuss below in detail, we find that the case 
of a H adatom on graphene is half a way between these two 
limiting cases. We show that the H impurity creates a defect 
with a complex magnetic structure where a spin S = \ is 
localized in a linear combination of the H and the neighboring 
C orbitals showing some of the properties of the vacancy-like 
defect. 

The paper is organized as follows: In Section [TT] we re- 
visit the Density Functional Theory (DFT) for diluted H on 
graphene and use the results as a guide to estimate realistic 
effective parameters. In Section [III] we present an Anderson- 
Hubbard model and its Hartree-Fock solution is compared 
with DFT results. In Section [IV] we show that the results can 
be interpreted in terms of an effective model of a H-C's clus- 
ter embedded in an effective medium and use the Numerical 
Renormalization Group (NRG) to analyze the stability of the 
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FIG. 1 . (Color online) 3D view of a supercell containing 72C atoms 
(large spheres) and one H impurity (small sphere) on top of a C atom 
(Co) represented by a green (light gray) sphere. The second nearest- 
neighbors C atoms (C n with n = 1,2,3) are represented by red (dark 
gray) spheres. 



magnetic solution. Finally, a summary and the conclusions 
are presented in SectionlVl 



II. DFT RESULTS 

We consider a graphene layer with a low concentration of 
H atoms. We calculate the spin dependent electronic struc- 
ture of systems with a unit cell of N C atoms and one H 
atom where N ranges between 32 and 72. Our DFT calcu- 
lations were done with a plane wave basis as implemented 
in the VASP code 44 and the cutoff energy was set to 400eV. 
The core electrons were treated with a frozen core projector 
augmented wave (PAW) methodi 4 ^ 4 ^ The frozen core pseudo- 
potential for carbon used a core radius for the s partial waves 
of 1.20 a.u. and 1.5 a.u. for the p channel. All partial waves 
for H where treated with a cutoff radius of 1.10 a.u.. We use 
the PBE generalized gradient approximation to treat the ex- 
change and correlations! 4 ^ 4 ^ To correct for the dipole moment 
generated in the cell and to improve convergence with respect 
to the periodic cell size, monopole and dipole corrections were 
considered^ 

Figure [TJ shows the relaxed structure of a supercell of 72 C 
atoms containing one H adatom. After minimizing the total 
energy with respect to the coordinates of all atoms in the unit 
cell, we find a H induced distortion consisting of a puckering 
of the hybridized carbon atom. The distortion is due to the 
modification of the electronic structure of the C atom bonded 
with the H adatom that changes from an sp 2 configuration to 
an sp 3 -like configuration after the hybridization with the H 
orbital. In the following, we use the label Co for the C atom 
directly bonded to the H impurity. Co is represented by a green 
(light gray) sphere in Fig. [TJ The second nearest-neighbors 
C atoms are represented by red (dark gray) spheres and are 
labeled as C„ with n = 1, 2, 3. 

Partial densities of states (PDOS) are calculated using a 
projector scheme provided in the PAW implementation of 
VASP4.6 when LORBIT=l 1 . The PDOS corresponding to the 
H atom shows the existence of a spin-split resonant state close 
to the Fermi level, see Fig.[2J a )- This induces a total magnetic 
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FIG. 2. (Color online) Partial density of states (PDOS) for the un- 
doped case projected on (a) the H atom, and the p z -orbitals of (b) the 
nearest-neighbor carbon atom (Co), (c) one of the second nearest- 
neighbors carbon atoms (C„), and (d) for a C atom located far from 
the adatom. For each case case both the majority and minority spin 
PDOS are plotted. A clear spin-polarized resonant state near the 
Fermi level Ef is observed in Figs.[2t a ) and|2jc). 



moment of « 0.49 /is/cell localized around the impurity. The 
covalent bond with the H induces a spin-polarization in the 
surrounding C atoms. The spin-polarization is small in the 
nearest C atom (Co), as shown in Fig. [2fb) where there is only 
a small spin-polarized peak close to the Fermi level. The mag- 
netic moment is larger at the second nearest-neighbor C atoms 
(C„) as can be appreciated in the strong spin-polarization 
of the p z orbitals of these atoms, see Fig. |2fc). The spin- 
polarization induced by the H adatom on the graphene layer 
tends to decrease at large distances from the absorption point, 
and far from the H adatom the PDOS corresponding to the p z 
orbitals of the C atoms [Fig. |2jd)] shows a Dirac point similar 
to the behavior of the C atoms in pristine graphene. 

The H adatom induces a magnetic moment and charge re- 
distribution on the surrounding carbon atoms that can be ana- 
lyzed by monitoring the atomic charge and magnetization as a 
function of the distance between the H atom and the graphene 
plane (R). This is shown in Fig. [3] The vertical dashed line 
marks the equilibrium position. The charge of the H atom in- 
creases, and its magnetic moment decreases, as it approaches 
the graphene sheet as shown in Figs. |3ja) and|3jb). The op- 
posite occurs for the second nearest-neighbors C„ atoms (red 
short-dashed lines) in Figs. [2c) and [3d), i.e., the charge of 
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FIG. 3. (Color online) Charge and magnetization of the most relevant 
atoms as a function of the distance R between the H atom and the 
graphene sheet. The equilibrium position of the adatom is indicated 
by a vertical dashed line. Notice how the charge (magnetization) of 
the H atom increases (decreases) as it approaches the graphene sheet. 
The opposite occurs for the C„ atoms (red short-dashed lines). Both 
the charge and the magnetization of the Co atom have a minimum 
close to the equilibrium position. 



the C n atoms decreases while the magnetization increases as 
the H atom approaches the graphene sheet. A different behav- 
ior can be observed at the absorption site, the Co atom. In this 
case, the charge shows a minimum and the absolute value of 
the magnetization a maximum close to the equilibrium posi- 
tion. 

Far from the absorption site, the perturbation induced by the 
H atom is smaller and alternates its magnitude between sub- 
lattices — it is larger on the sub-lattice corresponding to the C„ 
atoms. For example, the charge of the third nearest-neighbors 
[C nn atoms in Fig. EJc)] is almost independent of the position 
of the adatom, whereas the magnetization of these atoms only 
changes slightly compared with the changes in the magnetiza- 
tion of the C„ atoms [compare the red short-dashed and blue 
dashed lines in Fig. [2d)]. This shows that the H atom gen- 
erates a charge redistribution mostly on the same sub-lattice 
of the absortion site and predominantly at the neighboring 
C atoms, those that were highlighted with different colors in 
Fig. Q] As we will see in the next sections, the behavior of 
the system can be described by a simple model that takes into 
account the e-e correlations only in these atoms and connects 
them with the rest of the graphene lattice described with a 
simpler approximation. 

As was found recently for the case of fluorine adatoms, 50 
the nature of the chemical bonding of adatoms on graphene 
can change strongly with electron and hole doping. When 
graphene is electron-doped, the Co atom bonded to the F atom 
retracts back to the graphene plane and for high doping its 
electronic structure corresponds to nearly a pure sp 2 config- 
uration, see Ref. rf50ll for details. The situation is different 
for the electron doping of graphene with H impurities. To 
simulate electron doping we add one electron per unit cell. 
The extra charge is compensated with a uniform charge back- 
ground. The results obtained are shown in Fig. [4] The most 
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FIG. 4. (color online) Same as Fig. |2]but doped with an extra electron 
in the unit cell. Notice the absence of spin-polarization in the system 
after electron doping. 



important difference between Figs. [2] and [4] is tne absence of 
spin-polarization in the system after electron doping. As we 
can see, a similar electronic structure with well defined peaks 
at the Fermi level is observed but neither the H atom, nor the 
Co or C„ atoms show spin-polarization — see Figs. H|a)-|4tc). 
This dependence between gate doping and magnetic moment 
highlights the delicate interplay between electron correlations 
and localization in graphene with chemisorbed adatoms. In 
the following, we use our DFT results as a guide to formulate 
a theory that goes beyond the mean field level. 



III. ANDERSON-HUBBARD MODEL 

In order to interpret the DFT results and give a more phys- 
ical picture of the magnetic structure of the defect we use a 
simple Anderson-Hubbard (AH) model given by 

H = H graph + %imp + "Hhyb , (1) 

where 'Hgraph describes the 7r-bands of the graphene sheet, 
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here c\ a creates an electron with spin a at site i of the 
graphene lattice, the first sum runs over nearest neighbors, 
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FIG. 5. (color online) Comparison between the DFT PDOS (dashed 
line) projected onto the p z orbitals and the Hartree-Fock (solid line) 
results obtained with the full AH model. In both cases we use a 72C 
atoms unit cell. The right bottom panel shows the spatial distribution 
of the magnetization on the C atoms obtained with the AH model for 
a unit cell of 288 C atoms where finite size effects are negligible. 




FIG. 6. Schematic representation of the MAHM. The H impurity, 
the Co atom and the state |A) are treated exactly while the rest of 
the C atoms are treated as a non-interacting reservoir represented by 
the self-energy Ea(w). The capacitive coupling to the |±) states 
is taken into account only through a mean field approximation (see 
text). 



Uc is the Coulomb repulsion in the carbon atoms and n ia = 
cJ CT c i(T is the number operator of site i. The H impurity, which 
is bounded to the Co atom located at site i = 0, is described 
by 

Himp = y^gjf h\h a + U H ftj/^/ij/ij. (3) 

where creates an electron with spin a at the Is orbital of the 
H impurity with energy eh and intra-atomic Coulomb repul- 
sion Uh- The impurity-graphene interaction includes a one- 
body hybridization V and a distortion-induced shift of the Co 
carbon energy £o, 

Hhyb = £0 ™0 - V J2( C L h * + h Uo*) 
a 

with no = Cq^c ^ + cJ^c ^. In what follows we assume that 
the hopping matrix elements iy with i, j ^ are all equal 
to t = 2.8 eV while toj = tjo is reduced by the distortion. 
The energy of the hybridized carbon orbital is given by Eo — 
A 2 (e s — E p ), where e s ~ —8 eV and e p = are the energies of 
the s and p carbon orbitals, respectively, while A is a constant 
that parametrizes the deformation of the sp 2 bonding of Co- 15 
This expression interpolates between A — for the sp 2 and 
A = 2 for the sp 3, configurations. From the DFT results for 
the PDOS of a H far from the graphene surface we obtain that 
eh ~ -1.4* and U H ~ 2.8i while we set U c = t. Guided by 
the DFT results we take t j = to = 0.6t, Eo = —Q-7t and 
V = 2t. 

Figure [5] shows a comparison of the DFT results for the 
PDOS projected onto the p z orbitals with those obtained with 
the AH model within the Hartree-Fock (HF) approximation. 
In both cases we used a cell of 72 C atoms to facilitate the 
comparison. There is a good qualitative agreement between 
both methods, which indicates that the chosen effective AH 



parameters are adequate for capturing the relevant aspects of 
the problem — it should be emphasized that we do not intent 
to fit the parameters but find a reliable range of values for 
them. The overall qualitative agreement include some fea- 
tures related to the finite size of the cell as, for example, the 
appearance of gap-like and sharp peak structures. This is a 
consequence of the interference effects introduced by the pe- 
riodic array of impurities (all in the same sub-lattice). Within 
our TB model, such an effect can be eliminated, without much 
numerical effort, by increasing the size of the unit cell. Fig. 
|5jd) shows the spatial profile of the magnetization on a cell 
containing 288 C atoms, which shows a triangular symmetry 
characteristic of an isolated impurity. 



A. Minimal Anderson-Hubbard model 

The AH model presented above can be further reduced to 
a much simpler model that still captures the relevant aspects 
of the problem, allows to consider a single isolated impu- 
rity, greatly reduces the numerical work, and reproduces to 
an excellent accuracy the results of the full HF approach ob- 
tained with large unit cells. We first note that within the HF 
approximation, the presence of the H impurity generates a 
small charge redistribution mainly on its neighboring C atoms. 
Namely, in the Co atom and the three next nearest-neighbors, 
C„ with n = 1, 2, 3 (see Fig. [TJi. Therefore, it is sufficient to 
limit ourselves to consider a small cluster embedded in an ef- 
fective medium where the energies of the p z orbitals are fixed 
and the occupation numbers are solved self-consistently in the 
cluster. This approximation is numerically simpler as the self- 
consistent equations can be expressed in terms of integrals of 
analytical functions. Moreover, based on this approximation, 
the reduced Hamiltonian can be treated using more powerful 
numerical tools like exact diagonalization or the Numerical 
Renormalization Group (NRG) presented in the next section. 
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To illustrate the procedure, let us first consider the one-body 
part of %. Due to the hexagonal structure of the lattice, the Co 
atom only couples to the symmetric combination of the p z or- 
bitals of its nearest neighbor carbon atoms (C„). We denote 
that state as | A) and the corresponding fermionic operator as 

c Acr = En cJict/v^ where n — 1,2,3 labels the C„ atoms 
(we will also refer to this state as Ca)- The other two orthog- 
onal linear combinations of the C„ p 2 -orbitals, denoted by 
|±), are not directly coupled to Co- Furthermore, because of 
the symmetry of the hexagonal lattice, the states | A) and |±) 



are not coupled by the rest of the lattice either. As a result, the 
one-body terms of H can be separated into three decoupled 
parts — this is schematically shown in Fig. [6] 

Hence, for calculating the properties of the H impurity, the 
Co atom and Ca, it is sufficient to consider a reduced Hamil- 
tonian for the reduced system (see Fig. |6]l and include the rest 
of the lattice as an effective self-energy Sa(^) — the calcu- 
lation is presented in the appendix. More explicitly, the one- 
body Green function of the reduced system can be written as 
G(u>) = [uil - Uriu)]- 1 with 



n r {oj) 



(en 
V 





V 

so 








£A + ^aM 










-s H 










-XL 



(4) 



where we have included the possibility that the C„ atoms have 
a different energy than the rest of the C atoms of the graphene 
lattice (and included the |±) states for completeness). So far 
this is an exact procedure. The addition of the Coulomb inter- 
actions Uh and Uc only in Co can still be treated in a simi- 
lar way, provided we use the appropriate method, as the rest 
of the system remains a non-interacting fermionic bath and 
it can still be represented by Ea(w). This is no longer true 
when interactions are included in the rest of the C atoms. We 
will argue, however, that for a qualitative understanding it is 
sufficient, and important, to take into account the Uc inter- 
action only on Co and the C n (n = 1, 2, 3) atoms. We will 
refer to this model as the Minimal Anderson-Hubbard Model 
(MAHM). 

At the level of the HF approximation, the Coulomb re- 
pulsion u c Z)t=i,2,3 (*Wt ~ \) (™4 - \) shifts the enei "gy 

of the states |A) and |±) preserving the form of the effec- 
tive Hamiltonian H r (uj) for each spin projection and then 
the propagators. The spin dependent self-consistent Green 
function G a of the system is obtained from Eq. © with 
the self-consistent energies e a H — Eh + Un{{nHs) — §), 
£o = e + U c ((noa) - §), and 



£ A = £ A + U C 



'A/ 



(5) 



with a = -a. We have tested the validity of the MAHM 
by comparing the PDOS of the H, the Co, and the C n atoms 
with those obtained with the full model which includes the 
Uc interactions everywhere. The agreement between both ap- 
proaches is excellent, provided the unit cell used in the later 
case is large enough for the finite cell size effects to be neg- 
ligible. It is also worthy to emphasize that the MAHM leads 
to a magnetic structure similar to the one shown in Fig. |3d), 
which indicates that much of the observed antiferromagnetic 
structure is related to Friedel oscillations. 

It is interesting to compare results of the MAHM with those 
of DFT For that purpose, we plot in Fig. |7Ja) the magnetiza- 



tion of the H, Co and C„ atoms as a function of the hybridiza- 
tion V for two different values of eh ■ These results should be 
compared with those in Fig. [3] We clearly see that the MAHM 
is able to capture the most relevant features of the DFT results. 
Namely, that the magnetic state of the impurity is somewhere 
in between a pure adatom state (V <C t), where the magneti- 
zation is mainly localized at the H atom, and a vacancy state 
(V — > oo), where a substantial amount of magnetization has 
been transferred to the C„ atoms (mainly dominated by the 
|A) state). For a quantitative comparison of the results we 
have to consider the fact that in the DFT approach the lat- 
tice is relaxed for each H-graphene distance and consequently 
other parameters, like £q and to, also depend on the distance. 
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FIG. 7. (color online) Site magnetization as a function of the hy- 
bridization V between the H and Co atoms obtained within the 
MAHM for two different values of eh: eh = 1.294 (solid lines) 
and eh ~ 0-7t (dashed lines) and for Uc ~ t (a) and Uc = (b). 
Notice that the simple model is able to capture the general trend of 
the spatial distribution of the magnetic structure (compare with Fig. 
©. The presence of Uc not only induces some magnetization of the 
C„ atoms for large V but can also lead to a re-entry behavior of the 
magnetization. 
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Figure |7tb) shows the same parameters but in the absence of 
the e-e interaction on the C atoms. Clearly, the Uc interaction 
plays an important role in the case of large V (vacancy-like 
state), being responsible of the re-entry behavior observed in 
Fig. Ha). 

Before discussing the effect of the interactions beyond the 
HF approximation, we note that for lu — >• the self-energy 
Sa(w) au> In \d\ — i-\/3 \u>\, and the spectral density of the 
| A) state presents a divergence at the Dirac point, —Im(g^.) = 
= — Im(l/[a> — £a (<*>)])■ This is precisely the vacancy state 
(projected onto the |A) state) that has been extensively dis- 
cussed in the literatures^ This singular DOS is what makes 
the | A) state unstable against the formation of a localized 
magnetic moment when the interactions are included. Con- 
versely, theS± (w) self-energies diverge at the Dirac point and 
the spectral densities of the |±) states show a pseudo-gap at 
low energies (see the appendix). In view of this, one can ex- 
pect that the main role of the Uc interaction at low energies 
will manifest through the weakest coupled state |A). There- 
fore, we only keep the Uc/3 (n A ^ — h) (^AJ. — |) P af t of 



the interaction between the C„ (n = 1,2,3) atoms in what 
follows and neglect the rest. 

Under these assumptions the reduced Hamiltonian 
(MAHM) describes a correlated three-site cluster, given by 
the H, the Co and the | A) orbitals, embedded in an effective 
medium with a pseudo-gap described by Ea(^) — the mean 
field charge interaction between the | A) and the |±) orbitals 
that tends to keep the charge neutrality of the |A) state can 
be included through an effective site energy e a . In the next 
section we analyze the properties of this model. 



IV. EXACT DIAGONALIZATION AND NRG RESULTS 



The reduced Hamiltonian has the form 



n R = n cluster + ^band + rl 



hyb ; 



(6) 



where 



J 



"Hciustor = e H n H + U H h].h^h i + e n 



u c ( "or - 2 



e A n A + — I riAt - 2 



C Aa C 0a + C 0o C Act) 



(7) 



and "Hband and H^yb describe a band with a pseudo-gap at the 
Dirac point and the coupling to the | A) state 

riband = Zvflafva (8) 



HHyb = £ U(flc A<7 + c\j va ) , (9) 

that can be rewritten as 

n hyh = V2t Y^Ulc^ + cij 0a ). (io) 

where = t v f vcr , creates an electron on a sym- 

metric combination of H's third-nearest-neighbor C atoms. 

It is instructive to start the analysis by looking at the many- 
body states of the isolated cluster. For undopped graphene 
(/i = 0), we take e A = 0, Uh ~ 3£ and Uc = t as in the pre- 
vious sections but take £o = in what follows to simplify the 
analysis (taking £o 7^ introduces a electron-hole asymmetry 
which is not relevant at this point). 

We diagonalize ^cluster in the different charge sectors for 
different values of eh and V. Figure [8] shows the regions of 
stability of the different charge states in the [eh, V] plane. For 
Eh ~ — Uh /2 the ground state is always a 3-particle state 
with spin S = 1/2. The region of stability of the magnetic 
states as a function of eh shows a narrowing for V ~ 2t. For 




e H I u h] -ln(V/t) 

FIG. 8. (Color online) Color maps showing the magnetization of 
the H, C o and Ca atoms as a function of the orbital energy of the 
H atom, en and the adatom hybridization V as obtained from the 
exact diagonalization of ^cluster- The right lower panel shows the 
magnetization for eh = —Uh/2 (note that here the magnetization 
of C n is plotted). 



V < It, the magnetic moment is localized mainly at the H 
orbital while for V > 2t it is transferred to the | A) state. Note 
that, as the hybridization increases, the spin is transferred di- 
rectly from the H atom to the | A) orbitals, in agreement with 
the DFT and HF results. This is more clearly seen in Fig. |8jd) 
where the magnetization of the H, Co and C n atoms is plotted 
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as a function of ln(V/t) for ejj = ~Uh/2 (this is equivalent 
to plot it as a function of the H-Co distance if an exponential 
dependence of V is assumed). In agreement with the previous 
results, the magnetic moment is transferred to the C orbitals 
as the hybridization increases. 

From these exact results we can now calculate the Kondo 
coupling constant J. This is done in the standard way, elimi- 
nating high energy states of the cluster through a Schrieffer- 
Wolff transformation — to get 



J = 2t z 



I'CT 



(ft |4»HCA.| ^) + {f |CA»HcLl 4) 

E v — Ef 

(11) 

where \a) is the degenerate ground state of the cluster with 
energy K§ = E^, and the \v) are excited states. Note that in 
the absence of electron-hole symmetry in the cluster there will 
also be a local potential scattering for the conduction electrons 
due to the cluster. 

Figure[9|a) shows the J-coupling in the [eh, V] plane while 
Fig. |9|b) shows its dependence on V for different values of 
Br- For small V <C t, the magnetic moment, which is mainly 
localized on the H atom, is weakly coupled to the graphene 
sheet and the Kondo coupling is small. For a fixed V, the 
Kondo coupling increases as eh approaches the charge degen- 
eracy lines where perturbation theory fails and so J diverges 
due to vanishing denominators in Eq. ( fTTl i. For large V > 2t 
the bond between the H atom and the Co atom is strong 
enough to effectively decouple both atoms from the rest of the 
system. The magnetic moment is transferred to the | A) state 
and the Kondo coupling becomes V independent. In the large 
V limit it is simply given by J(V -> oo) = 8* 2 /^A = 24*. 

We complete the analysis by coupling the cluster to the 
rest of the graphene layer which provides an hybridization 
Ta ~ v3|w| at low energies. In the parameter region where 
the cluster is magnetic, we end-up with a pseudo-gap Kondo 
problemJ2iS^£ In this model, for an electron-hole symmetric 
impurity, the magnetic moment of the impurity (the cluster in 
the present case) remains unscreened even at zero tempera- 
ture. For the electron-hole asymmetric situation the magnetic 
moment can be screened only if the Kondo coupling is larger 
than a critical coupling J c of the order of the bandwidth3£i 
We calculated the stability of the magnetic moment in the 
cluster when it is coupled to the rest of the system. We used 
the NRG 63 ' 64 with a pseudo-gapped density of states and a 
Fermi velocity chosen to reproduce the effect of £a (u>). 

The solid line in Fig. |9ja) shows the NRG results for the 
stability region of the total magnetic moment on the cluster 
at zero temperature. We observe that the coupling to the rest 
of the system reduces the parameter area where the magnetic 
moment is stable at zero temperature (a similar effect is ob- 
served in the Hartree-Fock solution). For V < 2* the magnetic 
moment is mainly localized on the H atom and it is screened 
when the Kondo coupling reaches a critical value of ~ 10* 
(this unconventional screening is associated to a change of the 
cluster's charge 61 ). For larger values of V the magnetic mo- 
ment is localized on the |A) state. At finite but large (> 2t) 
values of V the critical coupling increases with increasing V. 
This is due to the fact that, in the present model, if V —> oo the 
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FIG. 9. (Color online) Kondo coupling constant J estimated from 
the exact diagonalization of the MAHM cluster, (a) as a function 
of eh and V and for the same parameters as in Fig. [8] (b) J as a 
function of V for three different values of eh' —Unfa (solid line), 
— 0.8£ (dashed line) and — 0.7t (short-dashed line). The case Uc = 
(dotted line) is included for comparison. 



magnetic impurity becomes electron-hole symmetric for any 
value of Eh and the magnetic moment remains unscreened for 
any value of J. 

It is worth mentioning that in the doped case the magneti- 
zation maps of Fig. [8] are modified^ 4 . For V > 2t, the mag- 
netic moment stability region narrows and shifts, being cen- 
tered around a line given by V 2 ~ a Eh + b. This shows that 
the impurity states is much more sensitive to doping in the 
vacancy-like regime than in the low V regime where the mag- 
netic moment is localized at the H atom. Therefore, a realistic 
H induced defect on graphene would be more sensitive to dop- 
ing than what one would expect from the usual Anderson-like 
model for an impurity on pristine graphene. A detailed study 
of this effect will be presented elsewhere. 



V. SUMMARY AND CONCLUSIONS 

We have combined a DFT description of diluted H impu- 
rities on graphene with tight-binding and effective models to 
describe the magnetic structure of a H induced defect. The 
DFT approach provides a realistic picture of the structural dis- 
tortions around the adsorbed H atom and predicts a magnetic 
moment localized in the neighborhood of the adatom. Ad- 
ditionally, it allows us to estimate parameters that are used 
to build an effective Anderson-Hubbard type model Hamil- 
tonian. The model is solved at the mean-field level and the 
results are compared with the full DFT band structure to test 
the quality of the mapping. The model Hamiltonian is then 
used to study a single H impurity adsorbed on graphene, a 
situation that can not be tackled with the current DFT based 
methods and that allows us to identify in detail the structure 
of the induced defect as well as its magnetic properties with- 
out the complications generated by the interactions between 
impurities. 

Within the single impurity Anderson-Hubbard Hamilto- 
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nian, the mean field approximation gives a magnetic solution 
for undoped graphene and a strong dependence of the impu- 
rity magnetic moment with doping — an effect with interesting 
implications for spintronics and for applications in magneto- 
transport devices. 

In order to treat the delicate balance between kinetic en- 
ergy and correlations at the defect including quantum fluctu- 
ations, we devised a minimal Anderson-Hubbard model that 
takes into account explicitly the electronic correlations at the 
impurity orbital as well as on the surrounding carbon atoms 
and replace the rest of the system by an effective medium. 

An analysis of the isolated cluster illustrates the structure 
of the magnetic moment. For —Uh < £h < and small hy- 
bridization V the spin is localized at the H orbital while for 
large V the spin is transferred to the carbon atoms forming a 
vacancy-like state. For intermediate values of the hybridiza- 
tion, that correspond to a realistic description of H, the stabil- 
ity region shows a neck and the magnetic moment is in a linear 
combination of the impurity and C orbitals. The effect of the 
rest of the host graphene is treated as an effective medium 
with a pseudo-gap using the NRG. We consider the case of 
undoped graphene where the Fermi energy lies at the Dirac 
point. As shown in Fig. [9] the region of stability of the mag- 
netic moment is narrowed as the cluster is coupled to the rest 
of the system. This behavior can be understood in terms of the 
known results for the Anderson impurity model in a system 
with a graphene-like pseudo-gap: in the case of electron-hole 
symmetry the spin is never screened, 54 while away from the 
electron-hole symmetry the spin can be screened at low tem- 
peratures if the Kondo coupling is larger than a critical value 
of the order of the bandwidth. Interestingly, within our model, 
in the large V regime (vacancy-like state), the electron-hole 
symmetry is recovered and the magnetic moment remains un- 
screened. 

While breaking the electron-hole symmetry will modify 
some of these results (a detailed study will be presented else- 
where), the main results of the present work are robust against 
it: (i) the spin is transferred to the carbon atoms as the hy- 
bridization increases, (ii) the Kondo coupling can reach quite 
large values, and (iii) for realistic values of the parameters 
(obtained from our DFT calculations) the H induced defect is 
half-way between the one corresponding to an adatom weakly 
coupled to pristine graphene and a carbon vacancy. 
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FIG. 10. Self-energy and unperturbed green functions of the C„ 
atoms, (a) and (c) correspond to the state [A) while (b) and (d) cor- 
respond to the |±) states. Note that Ea(w) presents a pseudo-gap 
close to the Dirac Point, and then <?a(w) shows a resonant-like be- 
havior (vacancy state). The opposite behavior is observed for E± (w) 
andg±(u). 



Appendix A: Non-interacting self -energy 

In Section [III] we solved the adatom problem using a sim- 
plified model where Coulomb interaction was assumed to be 
relevant only on the Co and C n carbon atoms. In that case, and 
if we are only interested on the properties of the reduced sys- 
tem formed by the adatom and the above mentioned C atoms, 
the presence of the rest of the graphene sheet can be taken into 
account through a self-energy contribution. In the following, 
we calculate this self-energy using the Dyson equation and 
taking advantage of the following: (i) the structure of the hon- 
eycomb lattice around a given atom has the same structure of 
a Bethe lattice (up to the second nearest-neighbors) which al- 
lows a simple disentanglement of the lattice Green function; 
(ii) the exact lattice Green function of a given site, Qi of the 
honeycomb lattice has an analytic closed form. 65 

Let us denote by = X)j c l/v^ tne fermionic opera- 
tor that creates an electron in an state that is a symmetric 
linear combination of the p z orbitals of the three nearest- 
neighbors atoms of Cq. The quantity of interest is the non- 
interacting self-energy of that state £^(w) due to the rest of 

the graphene (without Co). That is, if <7a(w) = ((c A ,cj^)) 
denotes the retarded Green function in the absence of the cou- 
pling with Co, then S^(w) = w — g~^{oS). The Dyson equa- 
tion, Q = g + gV Q, relates the unperturbed Green function 
g with the Green funtion Q in the presence of the perturba- 
tion V . Taking the hopping between the site "0" and its three 
nearest-neighbors as the perturbation, we can immediately ob- 
tain the following expression for the Green function of the site 
"0" (corresponding to the p z -orbital of Co), 



Go = 90 + 90 it Go #A 



(Al) 



9 



where g (oj) = (oj + iO+) _1 and 66 



2irt 2 V 2t 



21 



(A2) 



with 



S(x) = 8(V2x + 3- l)"5(V2a; + 3 + 3)"3 . (A3) 



and 



n , / K{k{x) 2 ) iflm(»lm(fc) 
yW_ i K(k(x) 2 ) + 2ip(x)K(l-k(x) 



<0 



\ fsT(A;(x) 2 ) + 2ip(a;)i ; i: (l-fc(x) 2 ) otherwise 

(A4) 

with K(k) the complete elliptic integral of the first kind, 

p(x) = sign(Imfx)) and k(x) = {2x + 3) 1/4 S(x)/2. Af- 
ter some straightforward algebra we finally get 



Sa(w) = lu 



?>i 2 



E (a;) 



(A5) 



where Eo(w) = w— C/ We notice that, in the low energy 
limit Cjcj-j ->■ 0), 

-Im(£ A H) ~ -3t 2 Im(£ (u;)) = y/3\u\ . (A6) 



This shows that the effective density of states of the part of 
the graphene sheet that couples to the symmetric state defined 
above also presents a pseudo-gap and goes linearly to zero 
[see Fig. [lOta)]. 

A similar analysis can be done for the green function and 
the self-energy corresponding to the |±) states, g±(ui) and 
S± (w), respectively. We then obtain 



3 1 



u) 



(A7) 



and 



E±(w) 



(A8) 



with Sa(^) = [w — Sa(w) — 3t 2 go( w )] 1 - This functions 
are plotted in Figs. [ToTb) andfTOtd). 

From these results, it is clear that, at low energies, the | A) 
state is only weakly coupled to the rest of the graphene sheet 
(excluding Co) while the opposite is true for the other two 
orthogonal states |±). This justifies our simplified model de- 
scribed in SectionUm 
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